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Abstract 

By using a reductive perturbation method, we derive from Maxwell-Bloch equations, a cubic 
generalized Kadomtsev-Petviashvili equation for ultrashort spatiotemporal optical pulse propaga- 
tion in cubic (Kerr-like) media, without the use of the slowly varying envelope approximation. We 
calculate the collapse threshold for the propagation of few-cycle spatiotemporal pulses described 
by the generic cubic generalized Kadomtsev-Petviashvili equation by a direct numerical method, 
and compare it to analytic results based on a rigorous virial theorem. Besides, typical evolution of 
the spectrum (integrated over the tranverse spatial coordinate) is given and a strongly asymmetric 
spectral broadening of ultrashort spatiotemporal pulses during collapse is evidenced. 

PACS numbers: 42.65.Tg, 42.65.Re, 05.45.Yv 
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I. INTRODUCTION 



The rapid progress over the past two decades in the area of ultrafast optics has led to 
the production, manipulation and control of pulses with durations down to a few optical 
cycles (see, e.g., the comprehensive review The theoretical and experimental studies 
of few-cycle pulses (FCPs) have opened the door to a series of applications in various fields 
such as light matter interaction, high-order harmonic generation, extreme nonlinear optics 



and attosecond physics ^j. On the theoretical arena three classes^ of main dynamical 



models for FCPs have been put forward: (i) the quantum approach 



fii) the refine- 



ments within the framework of the slowly varying envelope approximation (SVEA}of the 



nonlinear Schrodinger-type envelope equations [9 



-bj, and non-SVEA models (l3|-ll7j. In 



media with cubic (Kerr-type) optical nonlinearity the physics of (l + l)-dimensional FCPs 
can be adequately described beyond the SVEA by using different dynamical models, such 
as the modified Korteweg-de Vries (mKdV) 13 1, sine-Gordon (sG) fl^ . or mKdV-sG 
It is worthy to notice that the physics of the (l+l)-dimensional FCPs 

161]; this quite general model was also 



equations 

is well described by the generic mKdV-sG equation 
derived and studied in Refs. 



19 



20|. 



Another class of conceptually important optical problems for which the SVEA approach 
does not apply in the femtosecond regime includes multidimensional spatiotemporal optical 
solitons (alias "light bullets"), formed by the competing diffraction, dispersion, and quadratic 
2l| or cubic j22[ nonlinearity. Thus for the adequate description of multidimensional FCPs, 
a non-integrable generalized Kadomtsev-Petviashvili (KP) equation 23| (a two-dimensional 
version of the mKdV model) was introduced for (2+l)-dimensional few-optical-cycle spa- 



tiotemporal soliton propagation in cubic nonlinear media beyond the SVEA [2^, |25|. Re- 
cently, by using a powerful reductive perturbation technique 26|, or a multiscale analysis, 
a generic KP evolution equation governing the propagation of femtosecond spatiotemporal 
optical solitons in quadratic nonlinear media beyond the SVEA was also put forward 27 1. 

Transparency of the medium is obviously required by soliton propagation. It implies 
that all transition frequencies of the medium are far from the typical frequency of the wave. 
They can be either well above or well below the latter. As shown in (1+1) dimensions 

n 

[141 ]. a wave frequency much lower than the resonance frequency corresponds to a long 
wave approximation and to a mKdV model, while a wave frequency much higher than it 
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corresponds to a short wave approximation and a sG model. Later, a generic mKdV-sG 
model has been derived for the case of two transition lines, one well above, and the other 
one well below the wave frequency 16|, and is the most general non-SVEA model proposed 



to describe FCP propagation in Kerr media 



18j. It is expected to remain valid in the general 



case, where two sets of resonance lines are present instead of two single transitions. Our 
ultimate goal is to generalize the generic mKdV-sG model to (2+1) dimensions. However, 
both long- and short wave approximations for a simple two-level model have not yet been 
rigorously derived. It might be relevant to consider the full set of atomic levels, and it would 
be necessary in order to compute the nonlinear coefficients in a quantitative way. However, 
as a preliminary approach and for the sake of tractability, we need to consider a simplified 
model, but a natural question arises, why we consider a two-level model and not, e.g., a 
four-level one? As written above, it has been previously shown jl^ . [igI that if we consider 
only two separate transitions, one with resonance frequency below the optical range, and 
the second one above it, the FCP propagation is considerably affected. The most relevant 
refinement of the two-level model is thus a 2 x 2 level one. Our preliminary computations for 
the model of two sets of two levels show that in the long wave approximation considered in 
the present work, both the dispersive and the nonlinear coefficients merely sum up. In fact, 
we expect that the model equation obtained for the complete model with arbitrary number 
of levels, when all resonance lines are well above the wave frequency, has the same form as 
the one derived for the two-level model, only with modified values of nonlinear coefficients. 
However, this is not the case when the resonance line is well below the wave frequency, since 
the pulse evolution involves the difference of populations between the levels in a nonlinear 
way. In this situation, we expect a more complicated evolution equation for FCPs when more 
than two levels are involved. This general model, which might be, in our opinion, the most 
relevant correction to the two-level model, is left for further investigation. On the other hand, 
the two-level system has many interesting features by itself: It is in some sense equivalent to 
a classical oscillator, as can be shown from the computation of the nonlinear susceptibilities 
28l | , or from the derivation of the KdV model for FCPs propagation in quadratic nonlinear 



media 



291 ]. Moreover, estimations of the value of the nonlinear coefficient n2 drawn from 



the classical model have shown a very good agreement with experimental data, as, e.g., in 



the case of the model developed in Ref. 



30|. 



Thus there are three main issues which deserve to be investigated: the first one, which 
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constitutes the aim of the present manuscript, is the long- wave approximation leading to 
the cubic generalized Kadomtsev-Petviashvili (CGKP) equation, which yields collapse. The 
second and the third worth studying open problems are, respectively, the short wave ap- 
proximation and the investigation of the full model with both niKdV and sG terms. On 
this ground, we consider that the study of the CGKP model, including the derivation of the 
CGKP equation from a two-level model, is an important step in the adequate description of 
(2+l)-dimensional FCPs. 

The aim of this paper is to derive and to study a cubic generalized Kadomtsev-Petviashvili 
partial differential equation, which describes the dynamics of (2+l)-dimensional spatiotem- 
poral pulses in cubic nonlinear media beyond the SVEA model equations, starting from the 
Maxwell-Bloch equations for a set of two-level atoms. This paper is organized as follows. In 
Sec. II we derive the generic CGKP equation and we study by adequate numerical methods 
the propagation dynamics, the nonlinear difraction and, in the case of anomalous disper- 
sion, the collapse of ultrashort spatiotemporal pulses. We calculate in Sec. Ill the collapse 
threshold for the propagation of ultrashort spatiotemporal pulses described by the CGKP 
equation by both a direct numerical method and by an analytical method, which is based 
on a rigorous virial theorem. In Sec. IV, the evolution of the spectrum (integrated over the 
transverse coordinate) is given and a strongly asymmetric spectral broadening of ultrashort 
pulses during collapse is put into evidence. Section V presents our conclusions. 



II. THE CUBIC GENERALIZED KADOMTSEV-PETVIASHVILI EQUATION 
AND ITS NUMERICAL COMPUTATION 

For a Kerr (cubic) medium a CGKP equation can be derived from Maxwell-Bloch equa- 



tions using the powerful reductive perturbation method 



26| . As said in the introduction. 




we consider a set of two-level atoms with the Hamiltonian 

(1) 

where Q = Ub — Ua > is the frequency of the transition. The evolution of the electric field 
E is described by the wave equation 

{dl + d'^) E = ^^dl {E + AttP) , (2) 




where P is the polarization density. The hght propagation is coupled with the medium by 
means of a dipolar electric momentum 

/ n \ 

(3) 

directed along the same direction x as the electric field, according to 

H = Ho-i2E, (4) 

and the polarization density P along the x-direction is 

P = iVTr(p/x), (5) 

where is the volume density of atoms and p the density matrix. Since, as shown in jl^ , the 
relaxation can be neglected, the density-matrix evolution equation (Schrodinger equation) 
reduces to 

thdtp=[H,p]. (6) 

Transparency implies that the characteristic frequency of the considered radiation (in 
the optical range) strongly differs from the resonance frequency Q of the atoms. Here, as 
explained in the introduction, we assume that cu^ is much smaller than Q. This motivates 
the introduction of the slow variables 

e being a small parameter. The delayed time r involves propagation at some speed V to be 
determined. It is assumed to vary slowly in time according to the assumption Q. The 

pulse shape described by the variable r evolves even more slowly in time, the corresponding 
scale being that of variable (. The transverse spatial variable y has an intermediate scale as 
usual in KP-type expansions {2^. 



26| . To this aim 



Next we use the reductive perturbation method as developed in Ref. 
we expand the electric field E as power series of a small parameter e: 

E = €Ei + e^E2 + e^Es + ..., (8) 



as in the standard mKdV-type expansions 26|. A weak amplitude assumption is needed in 



order that the nonlinear effects arise at the same propagation distance scale as the dispersion 



does, which justifies that the expansion of E begins at order e. The polarization density P 
is expanded in the same way. 

The expansion [Eqs. ([7]) and ([HD] is then reported into the basic equations [Eqs. ([HE])], 
and solved order by order. The computation is very close to the (l+l)-dimensional case (see 
Ref. in what concerns nonlinearity and dispersion, while the treatment of dispersion 

and the dependency with respect to the transverse (spatial) variable rj is fully analogous to 
the case of quadratic nonlinear media, see Ref. j^^]. As a result we get the model equation 

d^drE, = AdtE, + {E,f + ^92^1, (9) 

which is a CGKP equation. The dispersion and nonlinear coefficients A and B in the 
above (2+l)-dimensional evolution equation are 



respectively, the refractive index being 



nchVt^ ' nclv'VL'^ ' ^ ^ 



The CGKP model equation ([9]) is generalized by expressing the coefficients A and B in terms 
of dispersion relation /c(a;) and third-order nonlinear susceptibility x^^'^ ■ They can indeed be 
written as 



, 1 d!^k 
A 



6 doj'^ 



Gvr /.IN 



(12) 
(13) 



w=0 



djjj I 



Also, the coefficient V in Eq. (|9]) is the group velocity, V — 

By means of an adequate linear rescaling, the CGKP equation ^ can be reduced to the 
normalized form 

[uz + aiU^UT + (T2Uttt)j. = uyy, (14) 

where ai = sgn{—B) and a2 = sgn(— A). In the frame of the Maxwell-Bloch equations, 
cTi = (T2 = — 1, hence the nonlinearity and dispersion yield temporal self-compression, but 
nonlinearity and diffraction tend to defocuse the FCP. 

The CGKP equation ([T3j) is solved by means of the fourth order Runge Kutta exponential 
time differencing (RK4ETD) scheme 311] • It involves one integration with respect to r. The 
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inverse derivative is computed by means of a Fourier transform, which imphes that the 
integration constant is fixed so that the mean value of the inverse derivative is zero, but also 
that the linear term is replaced with zero, i.e., the mean value of the function uyy is set to 
zero. For low frequencies, the coefficients of the RK4ETD scheme are computed by means 
of series expansions, to avoid catastrophic consequences of limited numerical accuracy. 
We consider in our simulations input data in the form 

u{T,Y,Z = 0) = A e-^'/P^-v-V?^ sin(a;T) . (15) 

Setting (Ti = (72 = — 1 as pertains to the Maxwell-Bloch system, we observe nonlinear 
diffraction (Fig. [1]). The nonlinear effect strongly increases the diffraction (compare Figs. [It 
and d). The conjugated effect of temporal self-compression and diffraction may lead to an 
intermediary stage, in which the pulse is very well localized temporally, and strongly widens 
spatially. This leads to a characteristic crescent shape (Fig. [lb). For the computations we 
used the parameters p = 4.0825, q = 2.8868, A = 4.8990, and u = 1. 

The CGKP equation (|9]) can be generalized using the expressions ( [T2|) and f [T3|) of the 
coefficients. For a medium with anomalous dispersion and focusing nonlinearity. A, B < 
and 0"! = (72 = +1- In this case spatiotemporal self- focusing occurs. 

Numerical resolution has been performed for u = —2.1909, p = 1.8257, q = 1.414, and 
for several values of the initial amplitude A, especially A = 3.80, 3.65, 3.57, and 3.50. Clear 
numerical evidence for collapse is found, see Figs. [2] and [31 Collapse occurs for the two 
highest values of the input spatiotemporal field amplitudes A (curves a and b in Fig. 2), 
and not for the two lowest ones (curves c and d in Fig. 2). Hence the occurence of some 
input amplitude threshold Ath is evidenced, and for the considered pulse shape, frequency, 
length and width, we get the numerical estimation 3.57 < Afh < 3.65. 

On the other hand, a rigorous mathematical analysis of the CGKP equation ([9]) has 
proved that wave collapse do occurs (for a comprehensive review of wave collapse in optics 



and plasma waves, see Ref. 



33| 



32|); also, for more details concerning CGKP equation see Refs. 



36| and the analysis performed in the next Section. 



7 



(a) 



(b) 
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FIG. 1: (Color online) Nonlinear diffraction of a FCP. a) Initial data {Z = 0). b) An intermediary 
stage, with a typical crescent shape {Z = 0.9594). c) Nonlinear diffraction at Z = 3.984. d) Linear 
diffraction at the same propagation distance for comparison. Input data is expression (|15p with 
p = 4.0825, q = 2.8868, u = 1, A = 4.8990 (a, b, c), and A = W'^ (d). 



III. CALCULATION OF THE COLLAPSE THRESHOLD 

In the following let us consider the generalized KP (GKP) equation in its normalized 
form 

{uz + U^Ut + uttt)t = '^YY- (16) 

The focusing CGKP equation Qj with ai = cr2 = +1 is obtained for the particular case 
p = 2. The GKP equation possesses the conserved Hamiltonian 

-// 

in which v = J u. Note that the existence of the conserved Hamiltonian H assumes that 
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1 2 1 2 

+ - 



1 

(p + l)(p + 2) 



u 



p+2 



dTdY 



(17) 




(a) (b) 

FIG. 2: (Color online) The collapsing pulse, a) Input; b) Last computed point. Input data is 
expression (fTSjl with to = —2.1909, p = 1.8257, q = 1.414, and amplitude A = 3.65 (curve b in Fig. 
[3] below). 
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FIG. 3: (Color online) Evolution of the maximal value of the electric field for a few values of the 
initial amplitude. The collapse occurs above some amplitude threshold. Input data is expression 
()15p with uj = —2.1909, p = 1.8257, q = 1.414, and several values of the amplitude A, namely 
A = 3.80 (a), 3.65 {h), 3.57 (c) 3.50 {d). 



both u and vy vanish at +oo and — oo, which implies that 

"+00 

udT = 0, 

i.e., that the mean value of the electric field is zero. 

Using a virial theorem, it has been proved that collapse occurs for p > 4 and H < 



FIG. 4: (Color online) Evolution of the width wy of the pulse, for the same initial amplitudes as 
in Fig. [3l Self-focusing occurs initially for all considered data. 



33|, |3J]. However, these conditions are sufficient 



the solitary wave is unstable for p > 4/3 
theorem was given in Ref. 



34 



3ut not necessary. Then, it was shown that 



35|, and a proof of collapse by using a virial 
36| for 4/3 < p < 4, which includes the particular case p = 2, 



which is relevant for the study of few-cycle optical pulses. The assumptions of the theorem 
involve the functional 

3p 



Q{u) 




2 2 
Vy + Ut ~ 



-u 



p+2 



dTdY. 



(18) 



2(p+l)(p + 2) 

Together with regularity conditions and conditions involving the solitary wave solution, 
the main assumption of the blow-up theorem is that the functional Q{u) < 0, in which 
u = u{T,Y,Z = 0) is the initial data. For the expression of u given by Eq. f[T^ . the 
functional Q{u) can be exactly computed using standard methods, and we get 



Q{u) 



128pq 



A^p^q^ ( 4e~ 



32g2 ( a;y 



l-e 



-F8p''V27re~ 



Im 



erf ( — 

\p 



tup 



dT 



(19) 



where Im is the imaginary part and erf is the error function. 

The factor which multiplies A^ in the wide bracket in Eq. ( fT9l) is always negative, and 
hence Q{u) is negative for A larger than some threshold value Ath-, with 



A 



th 



32q^ ( u; V 



Bp^v27r e" 



r°"im 

J —oo 



erf(f + T) 



dT 



p'^q^ (3 e-'^'p' - 4e' 



(20) 
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For the specific values of parameters used in our numerical computations, we find that 
Ath = 7.567. The threshold Ath — 3.6 found numerically is about half of this value Ath, 
found using the assumptions of the above mentioned virial theorem, which is consistent with 



the well-known fact that the conditions of the theorem of Ref. [36|], especially Q{u) < 0, 
yield only a sufficient condition. The quantity Ath will be referred to below as the 'analytical 
threshold for collapse' for the sake of simplicity. 

However, it is likely that this condition is optimal for an input which has a shape adapted 
to the collapse process. On the basis of this idea and of the numerical results, the discrepancy 
can be physically interpreted as follows. Figure S] presents the evolution of the width wy of 
the pulse against propagation distance Z, for a few initial values of input amplitude A (see 
Eq. f[T^ ). The width wy is numerically computed according to: 

2 _ / u^Y^dTdY 
Ju^dTdY ■ 

It is seen that self-focusing occurs in every case presented in Fig. 4. However, for initial 
amplitude below threshold At^, the self-focusing stops after a while and the collapse is 
inhibited. This feature is due to the dispersion (both linear and nonlinear), which tends 
to increase the temporal length of the pulse at the same time as it self-focuses. Below the 
threshold, the dispersion dominates and collapse is prevented, while above the threshold, 
self-focusing dominates and collapse occurs. It is worthy to mention that the arrest of 



24| 



collapse due to dispersion was already mentioned in Ref. 

This observation may justify qualitatively the discrepancy between the analytic and 
numerical thresholds for collapse found above. Between the two values for threshold 
[Ath — 3.6 ^ ^ ^ 7.6 ~ Ath), at the beginning of the process, the amplitude is not 
properly speaking sufficient to initiate collapse but, due to the shape of the pulse, a nonlin- 
ear lens effect induces a transverse self-focusing of the pulse, which increases the maximal 
pulse amplitude. At the same time, dispersion (linear and nonlinear) occurs, which tends 
to decrease the amplitude. If dispersion dominates, the growth of the amplitude stops and 
collapse does not occur. If, on the contrary, self-focusing dominates, the peak amplitude 
reaches a value which is sufficient to induce the collapse as such. Notice that in Fig. [3l the 
collapsing curves show two distinct parts, the first one (oscillating) corresponds rather to 
self-focusing and the second one corresponds rather to collapse stricto sensu. The value of 
the amplitude at the boundary between the two regions is close to the threshold value for 
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FIG. 5: (Color online) Evolution of the spectrum integrated over Y during collapse. The input 
data correspond to the case b in Fig. [3l 

collapse (Ath — 7.6) found from the mathematical condition Q{u) < 0, see Eqs. (fT9|) and 

m- 

IV. SPECTRAL BROADENING 

The evolution of the spectrum (integrated over Y) is shown in Figs. [5] and |6l At the 
beginning of the evolution process, the spectrum tends to be narrowed. Then the spectrum 
presents the oscilatory structure typical of the modulation instability (Fig. |6]d). However 
it is strongly asymmetric, in contrast with the case of the 'long' pulses described within 
the SVEA. Finally, a strong broadening of the spectrum is seen as the collapse itself occurs 
(Fig. Et). 

V. CONCLUSIONS 

In conclusion, we have introduced a model beyond the slowly varying envelope approxima- 
tion of the comonly used nonlinear Schrodinger-type evolution equations, for describing the 
propagation of (2+l)-dimensional spatiotemporal ultrashort optical solitons in Kerr (cubic) 
nonlinear media. Our approach is based on the Maxwell-Bloch equations for an ensemble of 
two level atoms and on the multiscale approach, and as a result of using the powerful reduc- 



tive perturbation method 26|], a generic cubic generalized Kadomtsev-Petviashvili partial 
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FIG. 6: (Color online) Evolution of the spectrum S (integrated over Y) during collapse, a) 
Initial data {Z = 0). b) At the limit between the 'self-focusing' and 'collapse stricto sensu^ stages, 
Z = 0.623. c) At the last computed time Z = 0.81445. The input data correspond to the case a in 
Fig. El 



differential evolution equation was put forward. Nonlinearly enhanced diffraction accom- 
panying temporal self-compression is observed. In the case of anomalous dispersion and 
focusing nonlinearity, collapse is evidenced. The collapse threshold for the propagation of 
ultrashort spatiotemporal pulses described by the cubic generalized Kadomtsev-Petviashvili 
equation was calculated numerically, and compared to the analytical results drawn from 
a theorem based on a virial method, which proves that collapse occurs. The discrepancy 
between the value of the collapse treshold obtained numerically and the corresponding an- 
alytical threshold for collapse is qualitatively explained. Moreover, the evolution of the 
spectrum (integrated over the transverse, spatial coordinate) is also given and a strongly 
asymmetric spectral broadening of ultrashort pulses during collapse is also put forward, in 
contrast to the case of long spatiotemporal pulses described within the slowly varying en- 
velope approximation. The present study is restricted to (2+1) dimensions, however, it can 
be extended to the (3+1) dimensions 37|] by incorporating into the generic model a second 
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transverse (spatial coordinate). 
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